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This study was a sensitivity analysis of a computer 
code (based upon CHAMPION/2/E/FIX), used at the Naval Post- 
graduate School mainly to model the flow within Solid Fuel 
Ramjets. The purpose of the study was to make the code 
"Simpler" and “more accurate" by reducing the required 
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SOurcCe term 

Stanton Number 

Temperature 

Axial velocity 
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Dimensionless distance of a point P from the solid 
boundary 


GREEK SYMBOLS 
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6 Incremental distance from the wall 
€ Turbulence dissipation rate; emissivity 
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eff 
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ref 


Combustion 
Convection 
Effective 
Flame 

Fuel 

Fuel grain 
Inlet 

Laminar 
Nitrogen 
Oxygen 
Near-wall node 
Products 
Reference 
Turbulent 
Referred to an independent property 
Wall 


Free stream 
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To EN ERODUC TION 


During the last fifteen years, there have been many 
advancements in the development of highly effective tacti- 
cal weapons for use at intermediate range and with high 
average flight speed. One of the simplest engines/motors 
for such use is the Solid Fuel Ram Jet (SFRJ). 

The SFRJ most often consists of a solid fuel grain which 
provides the walls for the combustion grain [Ref. lJ]. A 
sudden expansion at the air inlet of the combustor can be 
used to provide flame stabilization by inducing recircu- 
lation into the air flow (Fig. 1.1). Combustion can be sus- 
tained throughout the grain provided that the step height 
(h) is sufficiently large. However, the larger the step 
height, the larger the losses of the inlet stagnation pressure. 

Due to the sudden expansion, at least two distinct flow 
fields are generated within the fuel grain (Fig. 1.2). In 
the first field, the recirculation zone, the flow is highly 
turbulent and usually fuel rich. This hot gas region pro- 
vides the energy necessary to sustain the combustion proc- 
ess which must occur further downstream. Downstream of the 
flow reattachment, a turbulent boundary layer develops and 
the combustion is diffusion controlled. A diffusion flame 
emanates from the recirculation zone and remains within 


the developing turbulent boundary layer between the 
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Figure 1.2. Effects of Sudden Enlargements on a SFRJ 


fuel-rich zone near the wall and the oxygen-rich zone near 
the central core. 

There have been continuing research efforts concerning 
the behavior of the SFRJ at the Naval Postgraduate School. 
Both mathematical and experimental efforts have taken place 
to determine the effects of design as well as operational 


wemraoles of the system to be obtained. 
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Il. DESCRIPTION OF THE COMPUTER MODEL 


A. INTRODUCTION 

At the NPS there have been two basic computer codes 
used to model the flow within solid fuel ramjets. The 
first one was based on stream function (wy) and vorticity (w) 
[Ref. 3]. It predicted with reasonable accuracy the velocity 
and temperature distributions but did not predict pressure 
distributions accurately and was not easily extended to 
complex geometries or to three-dimensional flows. The second 
computer model is based on the "primitive" variables, pres- 
sure (p) and velocities (u,v) [Ref. 2]. It has also been 
used to predict the effects of fuel properties on the SFRJ 
performance and to evaluate the effects of different 


geometries as well as operating conditions. 


B. OVERVIEW OF THE PRIMITIVE VARIABLE PROGRAM 
1. Assumptions 

The CHAMPION/2/E/FIX computer program developed by 
Pun and Spalding [Ref. 2] has been used as the basis for the 
primitive variable model. The flow has been assumed to be 
steady (solid fuel grain dissipation is negligible), sub- 
sonic, recirculating, axisymmetric and to have constant 
specific heats. In order to calculate the effective viscosity, 


a modified Jones-Launder [Ref. 4] two-parameter turbulence 


16 





model has been used. It uses five empirical constants 
(Table II) and requires that two additional variables 
(turbulent kinetic energy (k) and turbulence dissipation 
rate (€)) be evaluated. The effective viscosity was calcu- 


lated using the formulas 


Mere ~ Hiam ~ Me. Peau? 
where 
2 
UW, = - o K /E (222) 
TABLE I 

k-e Turbulence Model Empirical Constants 

=) “> <p “x eff °c, eff 

43 ero 2 09 les0 iS 


The combustion was considered to be mixing-~limited 
(infinitely fast kinetics) with a simple one-step chemical 


reaction of the form: 


1 Kg fuel + i Kg, oxidizer ~*+ (1+i)kg products (2b oi) 


Four species were concerned: oxygen, nitrogen, fuel 


and products. Because of the infinitely fast kinetics 


ay, 
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assumed, no oxygen could exist at the fuel surface and at 
the surface was assumed to be isothermal. The turbulent 
Prandtl and Schmidt numbers were taken to be equal to unity 
and therefore the turbulent Lewis number and laminar Prandtl 


number were also taken to be equal to unity. 


C. GOVERNING EQUATIONS 
With the above assumptions, the governing equations for 


axisymmetrical flows can be written as follows: 


a 12° 3 (p 34 _12 7 2¢ _ 
7x PU) i a yp (Prve) rae Ox x eG ar? =f et (2.4) 
where » stands for the independent variables u, v, k, e«, 

he, Me, ~ MQ ,/4 and ny The first two terms are the "convec- 
tion terms", the next two the "diffusion terms". rs is the 
effective exchange coefficient for turbulent flow (Table I) 

and S, is the "source term" (Table II). The stagnation 


® 
enthalpy for unity Prandtl and Schmidt numbers is 


h = h+ (u*+v*)/2 + « (2.5) 

where for nonreacting flows 
he doe. T C27 6)) 
p ) 


and for reacting flows 


eg 


i 





= AH 
h = Ee a +- oe a Tor) Cer) 


Temperature has been calculated using Eqs. (2.5) 


and (2.7). Density has been calculated using the perfect 
gas law 
=) ee 
p RT (2.8) 


Since the quantities h (stagnation enthalpy), Me 
2 


and x = Me, Mo/ 1 


have identical governing differential equations and, in 
appropriate dimensionless form, identical boundary condi- 
tions, a representative equation had to be solved only once. 
The Meneionless enthalpy was selected to be 

ae 


H = = (2.9) 


in fg 


a I 
ry 


where he, is the stagnation enthalpy deep in the fuel grain 


and h. is the stagnation enthalpy of the air at the inlet 
dump plane. The corresponding dimensionless forms for the 
other quantities were: 

“N oer "No 


mh = _—42in 2 (2a) 
2 


=n 
"No in Nofg 


and 
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: x = cre 
2 Onn 


fg in 


ae Ul Ox 
(2.12) 


fu Ox 1 


1. Boundary Conditions and Solution Procedure 


Fixed boundary conditions (inlet velocity profile, 
et al.) were specified at the desired or experimentally 
determined values. Specified gradient boundary conditions 
were handled by setting the appropriate convection/diffusion 
coefficients to zero in the finite difference equation 
("breaking the link") and then entering the appropriate 
gradient through linearized "false" source terms. The. 
geometry as well as the appropriate boundary conditions are 
summarized in Figure (2.1). 

"Plug flow" was assumed at the air inlet, although 
this was not a computer program limitation. The turbulence 
kinetic energy was assumed to be uniform, with a value which 
Gorresponded to the approximate turbulence intensity of the 
inlet flow. Radial and axial gradients were set equal to 
zero on the centerline and at the exit, respectively. All 
the nonreacting solid boundaries were considered adiabatic 
with both u and v velocity components equal to zero (no-slip 


condition). 
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For simplicity, the boundary layer was assumed to be 
composed of a laminar sublayer which is then developed to 
turbulent without assuming a "buffer zone". For the 


laminar sublayer, 


+ - 
eI = ZS 
2 Yp ( ) 
where 
u ya OY 
a = __P_ and yr = ak 
p vt 70 p 


and for the rest of the layer 


+ 
£n (Ey_) 
ae = (eee (2.14) 
p K 
For E = 9.0 and «x = 0.4 (values found from experimental 
data) 
+ + 
u SOS a are ken 2k 
p te ( ) 


Solving the above two equations, it is found that 
transition from laminar to turbulent flow occurs at y. = I11.5. 
The value of Yn at each near-wall node is found from the 


relation: 





ss 
W 
(—) (2.16) 
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+ 
For eee 5 
Yp 


Milam“ 
1 = ele C27) 
W 6 
+ 
and for i mae 5 
- al/e 
iy eat <1, O KD (2.18) 


The value of ie was assumed to be uniform from the wall to 


ene hnear-wall grid point. So, 





ot: if oll 4 0 ile 6 (2.19) 
P = P Milam 


Then, the shear stress is calculated using the formula 


[Ref. 3] 

1/2 _ l/r 2 
ie Cy 0 Kp = p 65 = “— 
p 

ae el / 2 

1/2 EP et at 

ee nic eee) Zee 0) 
pb P lam 


where, as noted previously, 


+ 
&n(Ey_) 
us = —* 


Due to the steep gradients of the properties in 
turbulent flows near the solid boundaries, the source terms 


for kK and € at the near-wall nodes were expressed in terms 
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of the wall shear stress. Also, we provides the boundary 
condition for the u and v equations. 
The turbulence dissipation rate fs at a near-wall 


node p is given by the relation [Ref. 3] 


ee ae 4 ets /.2 
an = Cy KD /kK6 = = F2.430 (22 21.) 
where the length scale is presumed to be proportional to the 


incremental distance from the wall, 


R = at K 6 (2922) 

It was found that when using a sudden expansion 
geometry in reacting flows the near-wall dissipation rate 
feeto be increased on the step face te, = co ©/048) and the 
grid spacing adjacent to the grain had to be fine (Vp <. JAS) 
in order to obtain a temperature distribution in qualitative 
and quantitative agreement with experimental data [Ref. 1]. 
Equation (2.17) implies that the wall shear stress is calcu- 
lated assuming a linear velocity profile when y" Samal lees Ds 
Therefore, a near-wall grid point can lie within the laminar 
sublayer but the source terms for x and e€ imply that Woee/ Ham 
is much greater than unity. This precludes Yn from being 
Significantly less than 11.5. 

For reacting flows, the boundary conditions for the 


dimensionless properties (Eqs. (2.9) through (2.12)) were zero 
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at the inlet and unity "deep" in the fuel grain. These 
properties were considered to have zero gradients on non- 
reacting surfaces. 

The assumptions made for reacting flows (unity 
Prandtl and Schmidt numbers, simple chemical reaction, con- 
stant specific heat and stagnation enthalpy defined in Eq. 
(2.5)) result in a general boundary condition for all 
"conserved" properties (> .) on a surface which has mass 
transfer given by the relation [Ref. 2] 


| 2%. 
Mow hy aE) bw! @Po, bw 7 o6 £9? Mas $2) 


where ~ may represent "eh. or xX = M,,, - Mo/i- 
A mass transfer conductance (g) is often defined 
from the, relation 


IO, 


a) 


od or ° bw ~ SN 7 Te Bw oe 


where ¢. is defined as the free-stream value of the property. 
Substituting Eq. (2.24) in Eq. (2.23), the mass 


transfer from the blowing wall is 


Sn G(on, = Pcbw?/ ‘ec _bw ~ eee (2.25) 


where 
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p ey 
Be) = zoe S20 (206) 
c,bw °c, £9 


represents the mass transfer parameter (or blowing parameter). 


The value of g can be approximated using 
g = (ou) | St Ze2 1) 


From Reynolds' analogy with unity Prandtl number, the 


Stanton number is 


e a 
an =. = =e Cee 
2 (ou°) 
p 
where Ce is the local friction coefficient. Thus, 
16 
all 
gS | (2 329)) 
p 


Using the Couette flow approximation for the boundary-layer 


behavior with mass transfer, 


g = g” &£n(l + BP)/BP (2.30) 
where 
* ° 
g = iim (9g) (278) 
BP+0 
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In this case, BP has been evaluated from the solution 


of the energy equation using 


hh 
BP = _ (Zo 2) 
Aw 7 ee 


The wall shear stress has been evaluated using 


Bamei2-17) or Eq. (2.20) and modified with Eq. (2.30) 


Gremeee ee Qn(l + BP)/BP (2753) 

where ce is the wall shear stress without wall mass addition. 
The mass transfer conductance (g) was evaluated 
msing Eq. (2.29). The wall heat flux (q") on all solid 


isothermal boundaries was evaluated using Reynolds' analogy 


-4"/ (h,, -h) = = (2.34) 
P 

Since the blowing rates were small for the solid fuel 
ramjet (typically BP < 2.0), k and € were evaluated using 
Eq. (2.4) and the terms presented in Table II, which incor- 
porate the empirical constants of Table I. 

Five variables, namely u, v, kx, ¢« and H (or) h were 
solved using Eq. (2.4) in finite difference form. The line- 
by-line iterative procedure employed upwind differencing 
and underrelaxation in order to promote convergence. On 


each radial line, the mass flow rate was calculated using the 
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local density. The "error" in mass flow (compared with the 
summation of "mass-in" at the upstream boundaries) was used 
to uniformly adjust the axial velocity over the entire line. 
Then, the pressure at all downstream locations was adjusted 
to approximately correct for the momentum imbalance created 
by the uniform axial velocity and a "pressure correction" 
equation was solved for each cell on the line. Finally, 
local cell velocity (axial and radial) and pressure were 


adjusted to satisfy cell-wise continuity. 


A 
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IIIf. ANALYSIS 


A. INTRODUCTION 

AS previously discussed, the purpose of this study was 
two-fold; first, to reduce the CPU time required to run the 
program, without considerably affecting the accuracy of 
results and second, to make the program output better fit the 
experimental measurements without considerably increasing the 
CPU time. The methods used in this optimization procedure 
were the following: 

a. Monitor the effects of the number of traverses in the 
tridiagonal matrix (TDMA). 

b. Monitor the effects of the number of iterations. 

c. Monitor the effects of the grid spacing close to the 
solid fuel grain wall, without changing the number of 
grid lines, in either the x or the r directions 
(constant number of nodes). 

d. Monitor the effects of changes in the number of grid 
lines in both the x and r directions (the number of 


nodes was changed). 


B. EFFECT OF TDMA TRAVERSES 

The number of TDMA traverses in the main program is denoted 
by the parameter NTDMA. If there is a predominant direction 
of flow, this is arranged to be in the x-direction. Then, 


TDMA traverses are made on constant x lines, in order to make 
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the adjustments necessary. Generally, the greater the number 
of TDMA traverses, the better the accuracy obtained but the 


more CPU time is consumed. 


C. EFFECT OF NUMBER OF ITERATIONS 

The number of iterations in the main program is denoted 
by the parameter LSWEEP. A high number of iterations (or 
sweeps) is needed in order to have small residual values in 
the variables used. Here again the same difficulties as with 
the number of TDMA traverses appeared. The larger the value 
of LSWEEP, the better the accuracy but the greater the CPU 
time required. It is expected that for a high number of 
iterations, the increase in accuracy will be small whereas 


the CPU time will be excessive. 


D. EFFECT OF GRID SPACING CLOSE TO THE SOLID FUEL GRAIN WALL 

In the original program, the grid lines spacing close to 
the physical boundaries was small. For small grid spacing 
close to the solid fuel grain wall, the non-dimensional length 
. at near wall grid points was much less than 11.5 and the 
flow was assumed to be laminar. The velocity profile in this 
case is linear, given by 

i 


uos«=«‘éS a 56 (3.1) 
e Yi am 





For dense packing of grid lines, the flow may be con- 


sidered as laminar for more than the first node close to the 
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Toe lieenis case, the velocity profile will consist of 
several straight lines before turbulence will be encountered, 
as shown on Fig. (3.1). This could result in inaccurate 
predictions. The spacing was changed such that the first 
constant-r grid line was spaced far enough from the wall 

in order to have yy > 11.5. In this case, the velocity pro- 
file will be a logarithmic curve. If this increased spacing 
resulted in good predictions, then the number of constant-r 
grid lines could be reduced, thus reducing the total CPU 


time. 


E. EFFECT OF CHANGES IN THE NUMBER OF GRID LINES 

Generally, the greater the number of grid lines in both r 
and x-directions, the better the accuracy expected but also 
the greater the storage as well as the CPU time required to 
execute the program. Since it is desirable (but not neces- 
Sary) to have a maximum ratio of 1.5 for distances between 
adjacent grid lines in the same direction and a maximum length 
to width ratio of 10.0 for any cell, it is preferable that 
both the number of constant-x and constant-r grid lines be 
changed (denoted by NX and NY respectively in the main program). 
In such a case, if the grid spacing close to the solid fuel 
grain wall could be increased--as discussed above--the number 
of lines in both directions could be reduced without affecting 
the accuracy of predictions. This would result in reduced 


requirements for storage and CPU time. 
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Figure 3.1. velocity Profile Close to Solid Fuel Grain 


wall 
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IV. RESULTS 


A. INTRODUCTION 

All the options discussed in the previous chapter have 
been tried in the modified CHAMPION/2/E/FIX computer program. 
The original program had been modified in the past at the 
Naval Postgraduate School to account for an aft-mixing 
chamber [Ref. 2], and for compressibility and radiation 
effects [Ref. 6]. This program calculates the flow in a 
SFRJ, whose dimensions are shown on Fig. (2.1). The input 
values used in the program are shown in Appendix B, which 
is a copy of the BLOCK DATA subroutine. 

The results of interest are plotted and discussed in 
the following sections of this chapter. The plots mainly 
refer to two x-constant grid lines after the first and second 


sudden enlargements at x = 202 mm and x = 312 mm, respectively. 


B. EFFECT OF NUMBER OF TRAVERSES 

The original program has been modified to run with three 
different values of TDMA traverses: NTDMA = 3, NTDMA = 5 
(originally used value) and NTDMA = 7. Nothing else had 
been changed. An attempt to run the program with NTDMA 
greater than 7 resulted in a repeated series of overflows 
after a relatively high number of iterations had taken place 


(LSWEEP > 20) and the program execution was terminated. 
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Pe eeaem@perature Distribution 


The temperature profile variations were negligible 
in this case, as shown on Figs. (4.1) and (4.2). Thus, by using 
a smaller number of traverses, the accuracy of temperature 
distribution prediction was not significantly affected. 
2. Effective Viscosity Distribution 
Here again, the effect of variations of effective 
viscosity, as shown on Figs. (4.3) and (4.4), were negligible. 
8S. Axial Velocity Distribution 
The axial velocity distribution after the first and 
second sudden enlargements are shown in Figs. (4.5) and 
(4.6), respectively. Once more, the number of TDMA traverses 


did not affect the accuracy. 


C. EFFECT OF NUMBER OF ITERATIONS 

The program was run for three different numbers of 
iterations; for LSWEEP = 100, LSWEEP = 200 (original value) 
and LSWEEP = 400, without any other changes. 

ies yy Distribution 

The values of Ye after the first sudden enlargement 

were essentially unchanged. Some differences appeared after 
the second sudden enlargement but these results did not seem 
to depend upon the number of iterations; the peak value of 
val for LSWEEP = 100 was between the peak values for LSWEEP = 200 
and LSWEEP = 400. In any case, the value of yy defies 
accurate prediction because of recirculating flow; the relatively 


high value of yy was an indication of high turbulence. Thus, 
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a change in the number of iterations did not significantly 
affect the a ais Sie Uc JON. 
2. demperaturesaastripbution 

The differences in the temperature distributions after 
the first sudden enlargement were negligible, as shown in 
Figure (4.8). The peak temperature, 2700°K, which was an 
indication of flame location, moved almost to the mid-radius 
point. Then, a relatively steep reduction of temperature 
occurred, and at the centerline its value was less than 
2000°K. 

After the second sudden enlargement, the temperature 
distribution differences were also negligible, as shown in 
Figure (4.9). The peak temperature now occurred at the center- 
line of the SFRJ, which meant that the enlargement effectively 
promoted the mixing of oxygen and unburned fuel. 

3. Effective Viscosity Distribution 

Figure (4.10) shows that the variations in effective 
viscosity distributions after the first sudden enlargement 
were almost negligible. All three curves coincided close 
to the solid fuel grain wall and then had a small separation 
as they iiesched the centerline. The same effect can be seen 
in Figure (4.11). These small differences in Wore close to 
the centerline did not really affect the calculations since, 
as will be seen later, the axial velocity distribution close 
to the centerline was amost uniform. As a result, the shear 


stress value was also low and was not affected by the variations 
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of Woeee From the above, it is seen that a reduction in 
the number of sweeps to one half of the originally specified 
value did not significantly affect the effective viscosity 
distribution. 
4, Axial Velocity Distribution 

As for temperature and effective viscosity distribu- 
tions, the variations in the distribution of axial velocity 
were also negligible for this range of iterations. Close to 
the solid fuel grain wall, the values of axial velocity were 
the same in all three cases. Close the centerline there was 
a separation, with higher values attained for higher numbers 
of iterations for both of the areas after the first and the 
second enlargement (shown on Figures (4.12) and (4.13), 
respectively). 

Here, it was evident that increasing the number of 
iterations increased the accuracy attained; although, 
for this range of sweeps, the variation in the predictions 


did not exceed 33%. 


D. EFFECT OF GRID SPACING CLOSE TO THE WALL 
il. Yn Distribution 
The grid spacing was increased close to the solid 
fuel grain wall, such that turbulent conditions would exist 
at the near wall grid points all along the combustion length. 
The values of vy have been plotted in Figure (4.14). Values 


higher than 11.5 indicate turbulent flow, whereas lower values 


indicate laminar flow. 
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As shown on Figure (4.14), increasing the grid 
spacing close to the wall resulted in the values of yy being 
generally increased along the region from the first to the 
second sudden expansion. After the second sudden expansion, 
the values of Yn that corresponded to increased grid spac- 
ing were generally lower, thus indicating less turbulence. 
However, the values were greater than 11.5 and the flow 
was handled as turbulent from the second expansion to the 
exit. 

fe Hemserature Distribution 

The temperature distribution was dramatically 
affected by the grid spacing close to the solid fuel grain 
wall (i.e., the distance from the first node to the wall). 
Figure (4.15) shows that the original program predicted a 
peak temperature rather far from the wall. On the other 
hand, the program with increased spacing from the wall pre- 
dicted peak temperatures close to the wall (the larger the 
grid spacing, the closer the distance of peak temperature to 
the wall). A peak temperature farther from the wall means 
a thicker boundary layer and, in turn, a more fuel rich 
mixture ratio. 

Figure (4.16) shows that the original program pre- 
dicted complete mixing in the aft mixing chamber. Peak 
temperature had been predicted to be at the centerline of 
the SFRJ, whereas with increased spacing the peak temperature 


occurred just behind the step wall. 
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3. Effective Viscosity Distribution 


As shown on Figure (4.17), the effective viscosity 
distribution close to the well was not greatly affected by 
the spacing of the grid lines. However, close to the center- 
line the effective viscosity value increased with increasing 
spacing, while the value of Hose close to the wall was slightly 
greater in the original program. 

The same happened after the second expansion, as shown 
on Figure (4.18), although the curves obtained were now 
somewhat different. Again, the effective viscosity from the 
Original program was greater close to the wall, while, as 
the centerline was approached, this value dropped again below 
the values predicted for increased grid spacing close to the 
ed 

4. Axial Velocity Distribution 

As shown on Figures (4.19) and (4.20), the teens 
program predicted much higher velocities than the programs 
run with increased grid spacing close to the wall. It was 
also noted that the velocity profile of the original on 
was typical of og inar flow, where the ones for increased 
Spacing were typical of turbulent flows (almost uniform 
velocity close to the core). 

In this case too, the same general effects in the 
distribution could be noted as in those for effective 


viscosity. The original program predicted higher values than 
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the programs with increased spacing and the velocity pro- 

files corresponded to laminar and turbulent flows, respectively. 
These large variations in the u-velocity distributions 

shown on Figures (4.19) and (4.20) were not surprising con- 

Sidering the variations in the temperature distributions. 

Since, for increased grid spacing, the temperature was lower 

in the core, the velocity was also lower than that obtained 


in the original program. 


E. EFFECT OF A WIDER GRID SPACING 

Since this case had similarities with the previous one, 
only one case was run; with increased grid spacing in both 
x and r directions. The y, vs. xX distribution has been 
plotted, as shown on Figure (4.21). It was found that this 
case was not of particular interest. Along the solid fuel 
grain wall Yn was Slightly increased. This in turn meant 
that predictions would not be as accurate as those obtained 
increasing the grid spacing close to the wall only. So, no 


further analysis took place for this case. 
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V. CONCLUSIONS 


From the Analysis of Chapter IV it can be stated that: 

a. The computer program was not sensitive to changes 
in the number of TDMA traverses, NTDMA. Thus, a lower 
value, such as NETDMA = 3 can be used without reducing the 
accuracy of predictions. 

b. As expected, the higher the number of sweeps, LSWEEP, 
the greater the accuracy. However, if the number of iterations 
exceeded 200, the increase in accuracy was negligible. For 
any number of sweeps between 100 and 200, the accuracy of 
prediction was found to be very satisfactory. 

c. The program was very sensitive in the grid spacing 
close to the solid fuel grain wall. For increased grid spacing 
in such a way as to have yy > 11.5 for most of the nodes 
adjacent to the solid fuel grain wall, the results of the 
program fit experimental data of Well [Ref. 5], as shown on 
Figure (5.1). Also, this increased spacing close to the 
wall permitted the reduction of the number of CElcdeEnes 
in the r and, consequently, the x-directions, with consider- 


able savings in preparation as well as CPU time. 
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VI. RECOMMENDATIONS 


In view of the results of Chapters IV and V, the follow- 
ing are recommended, regarding the execution of the modified 
CHAMPION/2/E/FIX program. 

a. The number of TDMA traverses should be reduced from 
five to three. 

b. The number of iterations may remain the same or be 
reduced to one half without significant reduction in the 
accuracy of predictions. The use of any number of iterations 
greater than 200 is not expected to promote the accuracy 
Significantly. 

c. While establishing the grid, the nodes which are 
closer to the solid fuel grain wall should be selected in 
such a way as to result in values greater than 11.5 (at 
least for most of them). This will increase the accuracy 
of predictions and will enable the investigator to use fewer 


grid lines in both the x and y directions. 
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APPENDIX A 


C/ASPI#5 JOB (3238,1085) , 'GENMIX',CLASS=C CHAOO 
C/*MAIN LINES=(12) CHAOO 
C/ EXEC FORTXCG, REGION=1024K CHAOO 
C/FORT.SYSPRIN1 DD DUMMY CHAOO 
C/FORT.SYSIN DD * CHAO0O 
C FFF HHT Hs Hd tH EHTS FFF FH FF HF FF CHAO0O 
C ANDREA MH SVYSHS AYTO TO FAKELLO TON XREIAZOMAI GIA 
TH THESIS CHAOO 
C F#HFFEF HEE EHF HE FH FF HF HF Ft HHH FE FHF CHA0O 
BLOCK DATA CHAOO 
COMMON CHAOO 
1/CASE51/RIN, IYIN, IYINP1,IXIN,IXINMI1 ,TAULW (42), 

XPUSLW (42) ,CTAULW, CHAOO 
1 CXPLW CHAOO 
1/CASE52/1XIN2,IYIN2,NXP1,NYR1,NXR1M1,NYRIMI1, 

NYRIM2, CHAOO 
1 CYPTW1,CTAUT1 ,CXPLW2,CTAUL2 ,HWALL2 CHAOO 
1/CASE1/UINLET,FLOWIN, RPIPE, XPIPE,FXSTEP,HINLET, 

HWALL CHAOO 
2/DNY/ DYG(42) ,CYV (42) , FV (42) ,FVNODE(42) ,R(42), 

RDYG (42) ,RDYV (42) CHAOO 
2,RSYG(42) ,RSYV(42) ,RV (42) ,RVCB(42) ,RVSQ(42), 

SYG (42) ,SYGTR (42) CHAOO 
2,SYV(42) ,SYVTRV (42) ,Y¥ (42) ,YV (42) CHAOO 
3/DNYONX/AE (42) ,AN(42) ,AP (42) ,AS (42) ,AW(42) ,C (42), 

D(42) ,DIFE (42) CHAOO 
3,DIFN (42) ,DUW(42) ,DIFW(42) ,DU(42) ,DV (42) , EMUE(42), 

EMUN (42) CHAOO 
3,EMUW (42) ,USTAR (42) ,VSTAR(42) CHAOO 
3,PHIOLD (42) ,RHCE(42) ,RHON (42) ,RHOW(42) ,S (42), 

SPRIME (42) CHAOO 
3, VQLUME (42) ,CONN (42) , CONE (42) , CONW (42) CHAOO 
4/DNX/DXG (47) ,DXU (47) , FAREA(47) , FU (47) CHAOO 
4,FUNODE (47) ,KOUNT (47) , RDXG (47) CHAO0O 
4,RDXU(47) ,RSXG(47) ,RSXU (47) , STORE (47) ,SXG(47), 

SXU (47) ,X(47) , XU (47) CHAO0O 

COMMON CHAOO 
5/DJPHI/ IEW(IC) , ILAST(10) ,IMON(10) , IXNY(10), 

IZERC (10) CHAO0O 
5,KSOLVE (11) ,KRS (10) ,RELAX(10) ,RSREF (10) ,RSSUM(10), 

TITLE (10) CHAOO 

COMMON CHAOO 
6/DO/CCHECK, DP, FLOWPC, FLOWST,FLOWUP, GREAT, IPLRS, 

IPREF,IPRINT CHAOO 
6, IWEEEP, 1X, LX1LNY, IX1LNY1,IX2NY2,IXMON,IXP1,IXPREF, 

IXU, LYMON, LYPREF CHAOO 
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6, JEMU,JH, JLAST, JP, IPP, JRHO,JU,JV,JVP1,KALLED, 


KINPRI , KEHOMU 


6, KTEST, LABPHI , LSWEEP, NSOLVE, NTRAV, KRAD 
6,1YF,IYFL,IYL,IYLM1,IYLP1,NTDMA,NUMCOL 
6,N) ,NXMA) ,NXM1,NXM2 ,NXYG, NXYP,NXYU, NXYV 
6,NYNYMA) ,NYM1,NYM2 ,PI,RSCHEK ,RSMAX, TINY 


6,IPFRS1,IPFRS2,KIN, 


KEX 


COMMON/PROP, EMUREF , PREFF (10) ,PRL(10) ,PRT(10), 


RHOREF 


COMMON/D2D1/ARSL (47,10) , RSLINE (47,10) 


COMMON/D2D2/U (1560) ,V(1560) , TKE(1600) , TED(1600), 


H(1600) ,PP(40), 


7P (1444) ,RHO(1600) ,EML (1600) ,RSF (40,40) ,T(1600) 


DIMENSION F(14200) 
COMMON 


9 /TURB/C1,O2,CD,SCRTCD,CD25,ECONST, CTAUTW,CYPTW, 


TAUTW (47) ,GENK (42) 


, 


CHA0O 
CHAO0O 
CHAOO 
CHAOO 
CHAOO 
CHAOO 


CHAOO 
CHAOO 


CHAOO 
CHAOO 
CHAOO 
CHAOO 


CHAOO 


9YPUSTW (47) , FACTKE, FACTED,JTKE, JTED,CAPPA,STANTW (47), 


HBULK (47), 


9RUN,CP,GC,HINR,HWALLR, PREFR, BARM, JT, HOHECK 
9/COMB/AWALL (42) , FLUX (42) , FLOWAL (42) ,BP(42), 


RR (42) ,TEMP1 (42), 


Sieur, Raat , HUD ,HTR, TREE, STC,HC,FN2ZIN,FOXIN, 


9FMCLFU, FMCLPR,FMCLOX, FMOLN2 


9 /CHEMCO/FFU (42) ,FOX(42) , FN2 (42), 
9FLFUIT,FLFLO?T, FLOXIT,FLOXCT,FLN2IT,FLN2)T, 
9PCEFU , PCECX, PCEN2, RRAVE, LINE 


REAL*§ TL(12) 
LOGICAL*1 LTG(3) 
COMMON 
D/OPLOT/TL 


D,KPLOT (11) ,NINT(11) ,PHI(40,40) ,CL(25) ,XP(47), 


SCALEW , SCALEH, LTG 


COMMON/RAD/TWALL , TFLAME, SIGMA, EPSIW, TRAD 


E,I RADS ,RAD1,FUEL (40) ,QRAD(40) ,QCON (40) , FMASF (40), 


VOLL (40) 


DIMENSION DIFS (42) ,EMUS (42) ,RHOS (42) ,CONS (42) 
EQUIVALENCE (CONS (2) , CONN (1) ) 


EQUIVALENCE (DIFS(2),DIFN(1)), 


(EMUS (2) , EMUN (1) ) 


GENERAL,FLOW PARAMETERS 


I SUF 


3.1415926/ 


c EQUIVALENCE (RHOS (2) ,RHON(1)) , (ARBAE,AREAW) 
EQUIVALNECE (RHOS(2) ,RHON (1) ) 
EQUEVALENCE (F(1),U(1)) 

DIMENSION A(42) ,B(42) 

EQUIVALENCE (A(1) ,AN(1)), (B(1),AS (1) ) 
ee ee 

DATA GREAT,TINY,PI/1.E30, 

C 
DATA RPIPE,XPIPE, UINLET, HINLET ,HWALL/ 

10.026924,0.4628,92.220,0.0,1.0/ 

C 

€ 


DATA HWALL2/1.0/ 
DATA IXIN, IYIN/2,7/ 
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CHAOO 
CHAOO 


hi 





e 


DATA IXIN,IYIN/2,7/ 
DATA KTEST/0/ 
DATA KTEST/1/ 


CHAPTER 2 -------------- GRID 


c 


ANAANA 


DATA NXMAX,NYMAX/40, 40/ 


DATA NX,NY/40,33/ 

DATA KRAD/2/ 

KRAD=2 MEANS CYLINDRICAL POLAR COORDINATES. 

IF ONE WISHES THE PROGRAM TO CALCULATE THE GRID, 
SET X(1) OR Y(1) NOT EQUAL TO 0.0 


Bae Bom , 201, -023,-.038, .055, .07333, .09166, 
HOS SRF li2332, 
Proc ome. Wee 7 oS 31, .20164, .21997, .2383,.25663, 
-27343,.28468 
peeves, 29 7 bo. 30055, .5026,.3043, .3053, 
Geeeveew noel, .o09G, .31185,.31522,.32022,.32772, 
.seeo2, £35892, 
DeScGeee, .20292,.00692, 41392, .42892, .44192, 
~4543,.4643, 
EOP sai Ore Ole, Oven 0.27 
Datos Y70e0,.000Gl, .00141, .00271,.00401, .00531, 
Oe oo, O07 3S. , 
AOU GGL a0 Onl, .012235, .01436, .01586, .01686, 
MU OL 79S, 01828, 
Bev homen 0 los o, Oho, .OLI00,.015910,.01925, 


CHU ey UZ, 02112, ,O2eie, .02315,.02415,202515, 
DU 2Z50sme U2ZGM6 0m 01260924 ,0.,0.,0.,0.,0.,0.,0.,0.,0./ 


DATA FXSTEp/1.5/ 


ee = —— = —— — VARIABLES 


C-- 


---- VARIABLES SOLVED BY TDMA 


DATA JU,JV,JTKE,JTED,JH,JPP/1,2,3,4,5,6/ 


~--- AUXILIARY VARIABLES 


DATA JP,JRHO,JEMU,JT,JLAST/7,8,9,10,10/ 


---- THE NUMBER OF DEPENDENT VARIABLES SOLVED 


DATA NSOLVE/5/ 
DATA KSOLVE/9*1,2*0/ 
DATA KRS/6*1,4*0/ 


DATA MEBMeWA .3,0.3)70-5,0.5,1.,1.,0.1.,01,0.5,1./ 
DATA TITLE/10*4HXXXX/ 


CHAPTER 4 -<------<------ PROPERTY DATA 


DATA RHOREF, EMUREF/1.403,3.034E-5/ 
BAPA PRL APRT/12*1., .9, =%7*1./ 
BAiAweRmnn ot) .0, led, 6*1.0/ 


DATA Cl1,C2,CD,CAPPA, ECONST/1.43,1.92,0.09,0.4,9.0/ 


EPSIW IS WALL EMISSIVITY,RAD1L IS EMP.CONST.FOR 
EMISSGAS/UNBURNED FUEL ON A LINE 
DATA SIGMA,EPSIW,TWALL,TFLAME, RAD1/ 


Gi/ 


CHAOO 
CHAOO 
CHAOO 
CHAO0O 
CHAOO 
CHAOO 
CHAO0O 
CHAO0O 
CHAOO 
CHAO 0 
CHAO 0 
CHAOO 
CHA0O0 


CHAO 0 


CHAOO 
CHAOO 


CHAO0O 


CHAOO 
CHAO0O 


CHAO 0 


CHAOO 
CHAOO 
CHAO L 
CHAO] 
CHAO] 
CHAOL 
CHAOL 
CHAOL 
CHAO 1 
CHAO] 
CHAOL 
CHAO] 
CHAOL 
CHAO] 
CHAOL 
CHAOL 
CHAOL 
CHAOL 
CHAOL 
CHAOL 
CHAQ] 
CHAO] 
CHAO] 
CHAOL 
CHAOL 
CHAOL 





KS". Oca 00 ,02.05,000,,3190.,3.0000/ 


@ IRAD=0,BP NOT MODIFIED FOR RADIATION 
Cc IR/D=l1 WITH RADIATION 
DATA IRAD/O/ 
e DATA IRAD/1/ 
CHAPTER 5 -------------- STARTING PREPARATIONS 


DATA IXPREF,IYPREF/2,2/ 

DATA IXMON,IYMON/22,24/ 

DATA KINPRI/0/ 

DATA F,DV/14200*0.0,42*0.0/ 
DATA FACTKE, FACTED/0.005,0.03/ 


CHAPTER 6 --------------- STEP CONTROL 

CHAPTER 7 -------------- BOUNDARY CONDITIONS 
DATA KIN,KEX/3,1/ 

CHAPTER 8 -------------- ADVANCE 


DATA NTDMA/5/ 
DATA RSCHEK/0.001/ 
mee toR 9 =-—-—---=----__ COMPLETE 


SUBROUTINE DENSIT CALCULATES VARIABLE DENSITY 
ON EACH LINE. 

RUN=UNIVERSAL GAS CONSTANT 

GP=C SUB. P 

PREFR IS THE REFERENCE PRES ADDED TO P(I) TO 
CALCULATE P ACTUAL 


CH Qrera are 


DATA FLOWAL/42*0.0/ 
DATA RUN,OP,GC/8314.0,1004.0,1.0/ 


Cc 
DATA HINP,HWALUR, PREFR, BARM/2.6570E06,4.826E05, 
3.152E05,28.964/ 
DATA HTR,FOXIN,FN2IN,STC,HC/2.300E)6,0.233, 
0.767,1.92,2.443E07/ 
DATA FMOLFU,FMCLPR,FMOLOX,FMOLN2/100.0,28.97, 
32 OM2 8:..0,/ 
DATA RHOF,EMUF/1199.,2.43E-05/ 
DATA TREF/600.0/ 
C 
C SET VARIABLE 'LINE' = IX ON THE LINE ON WHICH 
C CONSERVATION CALCULATIONS DESIRED. 
DATA LINE/29/ 
C 
eee ioe eo --—-—----—- - ADJUST 
Sonn == —— - PRINT 


DATA NUMCOL/10/ 
DATA IPLRS,IPRINT/200,200/ 
C ### CORRECTION ON BASIC PROGRAM ### 
DATA IPLRS,IPRINT/200,200/ 
CHAPTER 12 ----------- DECIDE 
DATA CCHECK/0.001/ 
DATA LSWEEP.200/ 
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CHAO 1 
CHAO] 
CHAO] 
CHAO] 
CHAO Ll 
CHAOL 
CHAO 1 
CHAO] 
CHAO] 
CHAOL 
CHAO] 
CHAOL 
CHAO 1 
CHAOL 
CHAOL 
CHAO] 
CHAOL1 
CHAOL 
CHAO] 


CHAO] 
CHAO] 
CHAOL 


CHAOL 
CHAOL 
CHAO 1 
CHAO] 
CHAOL 


CHAO] 
CHAO] 


CHAO] 
CHAO] 
CHAO] 
CHAO] 


CHAO] 
CHAO 1 
CHAO] 
CHAO] 
CHAOL 
CHAOL1 
CHAO 1 
CHAOL] 
CHAO] 
CHAOL 
CHAOL] 
CHAOL 





CHtFFFFFFHFFHHEE EE EEE HH EE ES FF tH FEE EE FF HH 
C DATA LSWEEP/2/ 
CH#FFFFFFFF FFF FFHHFFFFFFFFHHE FES SES SHH SE + HH FH HF 
EC 
END 
PROGRAM MAIN (OUTPUT , TAPE6=OUTPUT) 
FOR CHAMPION 2/E/FIX,CASES50, (HTS7628,C50) 
IDENT CASE51 
29.7 MT 


/JULY 1977/ 


A DETAILED DESCRIPTION OF THE PROGRAM CAN BE 
FOUND IN THE REPORT A GENERAL COMPUTER 
PROGRAM FOR TWO-DIMENSIONAL ELLIPTIC FLOWS 
BY W.M. PUN AND D.B. SPALDING, HTS/76/2, 
MECHANICAL ENGINEERING DEPARTMENT, IMPERIAL 
COLLEGE, LONDON 


ANMAANANANANAANAANNANAINAANAN 


CHAPTER 0 0 0 0 0 0 0 0 #0 DECLARATIONS 0 O QO 
e 

COMMON 

1/CASE51/RIN, IYIN, ILYINP1, IXIN, IXINm1, TAULW(42) , 
XPUSLW (42) ,CTAULW, 

1 cCXPLW 

1/CASE52/IXIN2,IYIN2,NXR1,NYR1,NXRIM1,NYRIM1, 
NYR1M2, 

1 CYPTW1,CTAUTI,CXPLW2 ,CTAUL2, HWALL2 

1/CASE1/UINLET,FLOWIN,RPIPE,XPIPE,FXSTEP,HINLET, 
HWALL 

2/DNY/ DYG(42) ,DYV(42) ,FV(42) , FVNODE (42) ,R(42), 
RDYG (42) , RDYV (42) 

2,RSYG(42) ,RSYV (42) ,RV(42) ,RVCB(42) ,RVSQ(42), 
SYG (42) ,SYGTR(42) 

2,SYV(42) ,SYVTRV(42) ,Y¥(42) ,YV(42) 

3 /DNYONX/AE (42) ,AN(42) ,AP(42) ,AS(42) ,AW(42) ,C(42), 
D(42) ,DIFE (42) 

3,DIFN(42) ,DUW(42) ,DIFW(42) ,DU(42) ,DV(42) ,EMUE(42), 
EMUN (42) 

3,EMUW (42) , USTAR(42) , VSTAR(42) 

3, PHIOLD (42) ,RHOE(42) ,RHON (42) ,RHOW(42) ,S(42), 
SPRIME (42) 

3, VOLUME (42) ,CONN (42) ,CONE(42) , CONW(42) 

4/DNX/DXG(47) ,DXU(47) ,FAREA(47) ,FU(47) 

4,FUNODE (47) ,KOUNT (47) ,RDXG(47) 

4,RDXU(47) ,RSXG(47) ,RSXU(47) ,STORE(47) ,SXG(47) , 
SXU(47) ,X(47) ,XU(47) 
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9/DJPHI/ IEW(10) , ILAST(10) , IMON(10) , EXNY(10), 


I ZERO (10) CHAO2 
5,KSOLVE (11) ,KRS (10) ,RELAX (10) ,RSREF (10) ,RSSUM(10), 

TITLE (10) CHAO 2 
COMMON CHAO2 
6/DO/CCHECK, DP, FLOWPC, FLOWST,FLOWUP, GREAT, IPLRS, 

IPREF, ITPRINT CHAO2 
6,ISWEEP, IX,IX1NY,IXINY1L, IX2NY2,IXMON, IXP1,IXPREF, 

IXU, IYMON, IYPREF CHAO 2 
6,JEMU,JH, JLAST,JP,JPP,JRH0,JU,JV,JVP1,KALLED, 

KINPRI , KRHOMU CHAO 2 
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